Dirac quasiparticles and spin-lattice relaxation in the mixed state 
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We present the results of quantum-mechanical calculations, using the singular gauge transforma- 
tion of Franz and Tesanovic, of the rate of planar Cu spin-lattice relaxation due to electron spin-flip 
scattering in the mixed state of high-Tc cuprate superconductors. The results show a non-monotonic 
temperature and frequency dependence that differs markedly from semiclassical Doppler-shifted re- 
sults and challenges the assertion that recent experimental observations of the rate of planar Cu 
and O spin-lattice relaxation in the mixed state of YBa2Cu307_i point to antiferromagnetic spin 
fluctuations as a better candidate for the elementary excitations of the superconducting state. 



PACS numbers: 74.20.-z, 74.25. Nf 



The nature of the low-lying excitations in the mixed 
state of d-wave superconductors remains an open ques- 
tion despite extensivajmaeaxch over the past several M^S-, 
Recent experimentalEJErBffl and theoretical studiesOtla'u 
have shown that the frequency dependence of the planar 
Cu and O spin-lattice relaxation rate (T^^) allows NMR 
to be used as a local probe of these low- lying excitations, 
offering a powerful test of theoretical treatments of the 
electronic structure in the mixed state. 

A generic feature of the NMR experiments is that 
T^^ shows a strong frequency dependence, increasing 
sharply with increasing resonance frequency, i.e. with de- 
creasing distance from the vortex cores. This frequency 
dependence has been qualitatively reproduced within a 
semiclassical model af Doppler-shifted Dirac quasiparti- 
cles by Wortis et alB (WBK).|JIawever, experiments on 
the O spins in Y|Ba2Cu307_J3'LlQ and on the Tl spins 
in Tl2Ba2Cu06-i-5El have shown that the relaxation rate 
has a non-monotonic dependence on resonance frequency, 
first decreasing as one moves away from the vortex cores 
and then increasing as one approaches the minimum reso- 
nance frequency, precisely where the Doppler shift goes to 
zero. This is in clear disagreement with the semiclassical 
picture where the Doppler shift of the linear dispersion of 
nodal quasiparticles is what gives rise to a non-zero local 
density of states at low temperatures. An observed Jign- 
monotonic dependence of {TiT)~^ on temperatuijeEloEl is 
also in disagreement with the semiclassical modelQ which 
predicts that {TiT)~^ will always increase with tempera- 
ture. Different theoretical predictions were made by Morr 
and WortiaJ and Morm who calculated T{'^ within a spin- 
fermion model in which the damping of antiferromagnetic 
spin fluctuations is determined by their coupling to pla- 
nar quasiparticles. At low temperatures they found that 
T^^ increases monotonically with resonance frequency, 
but they also found that as the temperature is increased, 
(TiT)~^ actually decreases at the high-frequency end of 
the curve and increases at the low-frequency end such 
that a minimum develops at a temperature-dependent 
crossover point. Morr and Wortis noted that this tem- 



perature and frequency dependence is different from that 
calculated within the Dirac quasiparticle model and pro- 
posed that these differences would allow NMR experi- 
ments to determine whether it is Dirac quasiparticles or 
strong spin fluctuations that govern the low-temperature 
behaviour in the superconducting state of the high-Tc 
cuprates. In a later paper, where Morr repeated the cal- 
culations of Ref. for the planar O spins, Morr further 
argued, based on comparisons with the experimental re- 
sults of Curro et ala, that the antiferromagnetic spin- 
fluctuation mechanism is the dominant contributor to the 
relaxation rate for both the O and Cu spins. 

In this paper we show that a quantum-mechanical 
treatment of Dirac quasiparticles leads to results that 
are qualitatively different from those of the semiclassi- 
cal treatment of Ref. ^. In particular, we flnd that Tf^ 
has a non-monotonic dependence on temperature and on 
resonance frequency. These results complicate the ex- 
perimental distinction between Dirac quasiparticles and 
antiferromagnetic spin fluctuations, and suggest that the 
experimentally observed relaxation rate outside of the 
vortex core can be largely explained (depending on the 
orientation of the vortex lattice) by the influence of the 
periodic vortex lattice on the Dirac quasiparticles. 

In what follows we confine ourselves to a single Cu- 
O plane, and assume that the electronic properties are 
entirely two-dimensional. Following WBK, we consider 
the (— i <-*■ — |) transition of the in-plane Cu atoms and 
write the relaxation rate as 

rr^(r) = i^(7e7n/»^)|(-i|/+(r)| - DP 



{f\A^S-{r) + J2 B±S-{r + 6,)\i) 5{E, -Ef)\, 

(1) 

where i and / are the initial and final many body states 
of the electronic system, r is the position of a copper 
atom and <5,y are the displacements to the four nearest- 
neighbor copper atoms. ( ) denotes a thermal average, 
7e7„ft2B_L w 3.06 X 10-19 erg, and A^jB^ « 0.8. We 
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FIG. 1: The square vortex lattice used in this paper, divided 
into two sublattices. Note that the x- and y-axes are chosen 
to lie along the node directions. 



neglect the nuclear Zeeman splitting in the i5 function, 
because it is so much smaller than that of the electron. 
The nuclear matrix element, with /+(r), simply gives a 
factor of 3. The spin operators 



5-(r) = *Ur)*T(r) 



the London model, neglecting variation in the gap mag- 
nitude and the local magnetic field, and treat the vortex 
cores as point objects. For large values of the Dirac cone 
anisotropy {ajj = vf/va, where vp is the Fermi veloc- 
ity at the node -and va is the slope of the gap at the 
node) Mel'nikovtll showed that the quasiparticle wave- 
functions are confined by the vortex lattice to lie along 
the nodal directions. Knapp et alx3 further showed that 
Mel'nikov's ID approximation could be used to calcu- 
late an energy spectrum for the quasiparticles in excel- 
lent qualitative-agreement with that calculated by Franz 
and Tesanoviolj. Knapp et al. also showed that the ap- 
proximation could be improved by adding small numbers 
of plane waves along the nodal direction in a "quasi-lD" 
approximation. We use this approximation here to ex- 
pand the quasiparticle wavefunctions around each node 
as a sum over plane waves. For example, at the node 
(/cf, 0) (in this paper we set the nodes to lie along cither 
the X- or j/-axis) we write 



(2) 



are expanded in terms of Bogoliubov quasiparticle oper- 
ators 




7A"ikT"/.«kT(r) - 7;^„ki<«ki(i-)], (3) 



where the sum is to be taken over the node /i, the band 
index n and the wavevector k, and N±_ is the number of 
planar Cu sites. 

For a d-wava-superconductor in the mixed state Franz 
and Tesanovict2l have used a bi-partite singular gauge 
transformation (see Fig. to show that, close to the 
nodes, the quasiparticle wavefunctions are best described 
by Bloch waves propagating through a periodic potential 
created by the combined action of the applied magnetic 
field and the periodic structure of supercurrents wind- 
ing around a lattice of vortices. Note that we assume 



(r) 



i(fca:— K)-r 



(4) 



where p_f is the Fermi momentum at the node, a is the 
spin and K are the reciprocal lattice vectors of the vortex 
lattice. Ky is the cutoff wave vector along the y (trans- 
verse) direction and K!^. is the cutoff wave vector along 
the X (node) direction and is taken to be smaller than 
Ky in the quasi-lD approximation. Note that for com- 
putational convenience we neglect the ky dependence, a 
good approximation at low energies. 

If we substitute the plane wave expansion for the quasi- 
particle wavefunctions into Eq. (|^) , and perform the ther- 
mal average we find that 



J 



Tr'(r)=^(7e7«fi')'^E E /(eM«kT)[l-/(V"'k'i)] 



^ /ink /i'n'k' 



X (5(e^„kT ^ ^M'n'k'i) 



A I +AB< cos 



KK' 



(5) 



Note that we have assumed that we are always close 
enough to the nodes that 

cos {a[{pf,^ - Pf,,^) + {k^ - K^) - (fc^ - K'^)]) 

^ cos {a[p^^-p^, x]), (6) 

where a ~ 3.855 A is the distance between adjacent 
copper atoms. In the calculations presented here, we 



have used the numbers in Chiao et aZJi3 for YBCO: 
Vp — 2.5 X 10^ cm/s, ao = 14. In the plane wave expan- 
sion of Eq. (H) we have used 45 reciprocal lattice vectors 
along the transverse direction and 11 along the node di- 
rection. 

The calculated spatial structure of-Tj^^ is shown in Fig. 
^. As in the semiclassical calculation^, the rate is highest 
near the vortex cores and falls with increasing distance 
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FIG. 2: (ri(r)T)"^ in a two-vortex unit cell 
applied field of 16 T. The peaks at the vortex positions have 
been cut off in the top panel to better show the spatial struc- 
ture. 



from the core. What is not seen in the semiclassical cal- 
culation is the complicated four-fold spatial structure of 
T^^. This structure is entirely a quantum effect which 
is due to the extension, at low energies, of the quasipar- 
ticle wavefunctions along the node directions and by the 
confinement of the quasiparticle wavefunctions by the pe- 
riodic potential of the vortex lattice. The degree of spa- 
tial anisotropy increases with the size of the Dirac cone 
anisotropy, ajj. Note that this structure is not due to 
variation in the gap magnitude at the vortex core, which 
we neglect here. The extension of quasiparticle states 
along the node directions is observed in other theptetical 
treatments of quasiparticles in the vortex latticdlj'E^lEj, 
and the analogue of these states for the case of a single 
vortex is discussed in Refs. p7|p^ , p^ . 

The position-dependent rate can be used to generate a 
frequency-dependent rate (shown in Fig. ^ by correlating 
the position in the vortex lattice with the local magnetic 
field, as calculated using the London model: 



B 



K 



iKi 



(7) 
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FIG. 3: (TiT)~^ as a function of local magnetic field, B (or 
resonance frequency, ui oc B) for H = 8 T (o), H = 16 T (o) 
and T = 9 K. Note the upturn in the rate at the minimum 
in the local field (frequency), and the large spread in T^^ . 
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where we have used A ~ 1600 A and ^ sa 16 A. As 



FIG. 4: In order to have an upturn in T^^ at low frequencies 
lines drawn from the vortices along the node directions must 
cross the local field minimum. 



noted already, the rate increases with increasing fre- 
quency, however there is also an upturn at the minimum 
frequency. This upturn is due to the intersection of the 
extended quasiparticle wavefunctions with the local field 
minimum (see Figs. |[ §). It is in this sense that the low- 
frequency dependence is sensitive to the orientation of 
the vortex lattice. In the case of a square vortex lattice, 
this effect will be seen if the lattice is oriented along the 
Cu-0 bonds. 

For applied magnetic fields below 6 T the 
tice is believed to be triangular or obliquec 
in which this effect might be seen for certain orienta- 
tions. Themis.|^sa,ph|Owever, strong experimentalnJ and 
theoretica]E§n3ua'E2lE2l evidence for a transition from a 
triangular to a square vortex lattice with increasing mag- 
netic field. While the theoretical results suggest that the 
square lattice should be oriented along the node direc- 
tions, thft-jrecent small-angle neutron scattering results 
for LSCCxJ see very clear evidence for a field-based tran- 
sition between 0.5 and 0.8 T to a square lattice that is 
oriented along the Cu-0 bonds. Due to the smaller gap 
in LSCO, Gilardi et al. suggest that a similar transition 
might be seen in YBCO for fields a factor of 10 higher. 
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FIG. 5: {TiT)~^ as a function of temperature at 16 T, aver- 
aged over a narrow frequency region near the core (solid line), 
at the minimum frequency (dashed line) and at the frequency 
where the rate goes to a minimum (dotted line). The inset 
shows a close-up of the slower rates. 



FIG. 6: (TiT)"^ averaged over the vortex unit cell as a 
function of applied magnetic field for T = 1 K (o), 9 K (A), 
and 19 K (o). The structure near 15 T is due to shifting 
of the lowest quasiparticle energy bands which is caused by 
increasing the magnetic field. 



on the order of 5 to 8 T. While the exact structure and 
orientation of the vortex lattice remain controversial, we 
would argue that there is good evidence for the type of 
lattice in which Dirac quasiparticles would contribute to 
a low-frequency upturn in T-^^. An excellent test of this 
theory would be a frequency-dependent measurement of 
T-^^ in a material known to have a vortex lattice which 
does not have the orientation shown in Fig. IJ. 

The temperature dependence of (TiT)^^ at three dif- 
ferent frequency positions in the unit cell is shown in Fig. 
^. At all positions, the relaxation rate starts from zero 
at r = and quickly rises to a peak near T — 2 K. 
Away from the vortex, (TiT)^^ then decreases until 
T ~ 10 — 20 K, where it starts to increase agam. This 
is in contrast with the spin-fluctuation resultsQEi where 
(TiT)~^ was found to increase with temperature at the 
low-frequency end and decrease at the high-frequency end 
of the spectrum. The temperature dependence due to 
Dirac quasiparticles can be understood from the quasi- 
particle density of states. Following from the work of 
Franz et al£3, Marinelli et allEl showed that the density 
of states starts from zero at zero energy, but increases 
very rapidly due to the existence of energy bands which 
are bent back toward the Fermi surface by the periodic 
potential of the vortex lattice. Low energy structure in 
the density of states created by the bending of theJDw- 
est energy bands was further shown by Knapp et aliS to 
give a specific heat which exhibits similar behaviour to 
the temperature dependence of {TiT)~^ shown here. We 
note that this temperature de»endencQ-|is similar to that 
seen experimentally for the (Ma and TH spin lattices and 
that, unlike the low frequency upturn in T'j""'^ ^ the tem- 



perature dependence is not dependent on the orientfition 
of the vortex lattice. 

From Figures and o one can see that increasing the 



applied magnetic field has two effects: (i) at larger fields, 
the local magnetic field becomes more uniform and the 
frequency range narrows, but the overall shape of T^^ 
stays the same (see Fig. |3|); (ii) the low- lying density 
of states scales with y/H due to the existence |a|j-4)car- 
nodes in the quasiparticle energy spectrumtj'EBEil, so 
that {TiTy^ ^ H (see Fig. |) in agreement with the 
experimental findings of Zheng et a/.EHl for ^■^Cu NMR in 
TlSrsCaCuaOg.s- 

In conclusion, we have shown that the quantum- 
mechanical treatment of Dirac quasiparticles is able to 
qualitatively account for non-monotonic structure in the 
frequency and temperature dependence of (TiT)~^ in the 
mixed state of high-Tc cuprate superconductors. The 
quantum-mechanical results differ markedly from the 
semiclassical results of Ref. ^ and show some similari- 
ties with the spin- fluctuation picture of Refs. These 
similarities complicate the experimental distinction be- 
tween Dirac quasiparticles and snin fluctuations and call 
into question the recent assertion^ that antiferromagnetic 
spin fluctuations are the dominant mechanism for spin- 
lattice relaxation in the mixed state. 
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